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We analyze how the measurements of the low-energy spectrum of cosmic antiprotons can provide 
information on relic neutralinos. The analysis is focused on the light neutralinos which emerge 
in supersymmetric schemes where gaugino-mass unification is not assumed. We determine which 
ranges of the astrophysical parameters already imply stringent constraints on the supersymmetric 
configurations and those ranges which make the antiproton flux sensitive to the primary component 
generated by the neutralino self-annihilation. Our results are derived from some general properties 
of the antiproton flux proved to be valid for a generic cold dark matter candidate. 
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I. INTRODUCTION 

In supersymmetric schemes, where gaugino-mass uni- 
fication is not assumed, the lower bound on the neu- 
tralino mass is determined by the upper limit on the 
contribution of cold dark matter (CDM) to the cosmo- 
logical density parameter, (^CDM)max^ 2 - Using the value 
(^CDM)maxfr = 0.13, derived from results of Refs. 
one obtains the lower bound m x > 7 GeV 0, • This 
is at variance with the more commonly employed lower 
limit of about 50 GeV, derived from the experimental 
(LEP) lower bound on the chargino mass. The rich phe- 
nomenology related to possible light neutralinos, with 
masses within the range 7 GeV < m x < 50 GeV, has 
been discussed in Refs. [H 0, H, @ ■ 

In particular, in Ref. || it was scrutinized the capa- 
bility of indirect measurements of WIMPs (Weakly In- 
teracting Massive Particles) to detect light neutralinos 
and it was concluded that measurements of cosmic an- 
tiprotons represent the most promising mean of indirect 
exploration. The analysis of the antiproton signal was 
based on a study of the antiproton flux at one single 
representative value of the antiproton kinetic energy (Tp 
= 0.23 GeV). It was proved that, for a wide range of 
the astrophysical parameters, the antiproton signal due 
to neutralino pair-annihilation is within the level of de- 
tectability at small values of m x . 

In the present paper we extend the previous investi- 
gation of the antiproton signal by analyzing the detailed 
features of the expected theoretical spectra, and discuss 
how our results can be employed to determine the pres- 



ence of an actual primary signal or, at least, to derive 
significant constraints on supersymmetric configurations 
at small m x . Our analysis is carried out in the perspec- 
tive of a significant breakthrough in the determination of 
some relevant astrophysical parameters and in view of a 
sizeable improvement in the measurement of the cosmic 
antiprotons spectrum, as expected in forthcoming exper- 
iments in space. 

This paper is organized as follows. In Sect. [H] we 
analyze some properties of primary cosmic antiproton 
fluxes due to self-annihilation of a generic candidate of 
cold dark matter. These features are derived in the stan- 
dard scheme usually employed to describe the decoupling 
of cold particles from the primordial plasma. An upper 
bound is obtained for the antiproton flux. Then, in Sect. 
IIIII a full evaluation of the antiproton spectrum is derived 
in the case of relic neutralinos. Results and perspectives 
are presented in Sect. II VI 

II. A FEW PROPERTIES OF THE 
SELF-ANNIHILATION CROSS-SECTION FOR A 
GENERIC WIMP 

Cosmic primary antiprotons can originate from the 
hadronization of quarks and gluons produced in WIMP 
self-annihilation processes (we consider here selfconju- 
gate WIMPs) 0, H, H EH EH Once antiprotons are 
produced in the dark halo, they diffuse and propagate 
throughout the Galaxy. The propagated antiproton dif- 
ferential flux at a generic point of coordinates r, z in the 
Galactic rest frame (r is the radial distance from the 
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Galactic center in the Galactic plane and z is the verti- 
cal coordinate) is 



section times the relative velocity of the self-interacting 
particles 



with 
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holds (x is defined as x = m x /T, T being the tempera- 
ture) . 

For relic particles in the Galactic halo x ~ 10 6 , then, 
usually, a good approximation is: 



Notations are as follows: Vp is the antiproton velocity, 
(cann^)o is the average, over the Galactic velocity dis- 
tribution, of the WIMP annihilation cross-section multi- 
plied by the relative velocity. S^ Btl . is a function which 
takes into account all the effects of propagation in the 
Galaxy and includes p 2 (r,z), p(r,z) being the galactic 
dark matter distribution. Here we take: 
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where pi is the total local dark matter density set at 
the value of 0.3 GeV/cm 3 , a is a core parameter, a=3.5 
kpc and i? Q =8 kpc. In the following we do not include 
any dumpiness effect. £ represents the fractional lo- 
cal density p x of our generic WIMP as compared to pi, 
i.e. t = p x / Pi- By applying the usual rescaling procc- 
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From 



dure |12j , one has £ 

the analyses of Refs. |2| one derives that at 2a level 
the cosmologically interesting region for cold dark mat- 
ter is: 0.095 < 51cdm^ 2 < 0.13 (in what follows this 
will be denoted as the WMAP range for CDM abun- 
dance). Thus, for (r2cDM^ 2 )min we use here the value 

(f>CDMfr 2 )mm = 0.095. 

dNp/dTp is the differential antiproton spectrum per 
annihilation event: 
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where F denotes the different annihilation final states, 
BR(xx —* P + -X") the branching ratios and dNp/dTp 
stands for the antiproton energy spectra in the F channel. 

Throughout this paper we will be interested in the an- 
tiproton differential flux at Earth (r = Rq , z = 0) as 
a function of the antiproton kinetic energy; this flux will 
be simply denoted as $p(Tp). 

We turn now to a discussion about some general prop- 
erties of the quantities entering the flux factor T. To 
make the discussion more transparent, we consider a sce- 
nario where the standard expansion in S and P waves for 
the thermally averaged product of the annihilation cross- 
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We recall that < <7 ann u > enters also in the relic abun- 
dance: 
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3.3 x 10" 38 cm 2 
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where < a arm v > = x f {a ann v) int , (<7 ann t>) int being the 
integral of < cr ann v > from the present temperature up 
to the freeze-out temperature Tf\ Xf is defined as xj = 
m x /Tf and g i ,(xt) denotes the relativistic degrees of free- 
dom of the thermodynamic bath at Xf . Using the expan- 
sion of Eq. ijjjj, one obtains: (cr alm v) ~ d + l/(2x/)b. 
Since Xf ~ 20, also the P-wave contribution b has to 
be retained in this case. In the following, we however 
specifically assume that d > l/(2xf)\b\. 



A. Lower bound on < <7ann« >o 

The cosmological upper bound Q. x h 2 < (f2cDM^ 2 )max 
implies that (cr aim v) is limited from below. For a cold 
relic with a mass in the range 10 GeV < m x < 1 TeV, 

1/2 

one has g+(xf ~ 20) ~ 90, so that Xf/g*(xf) ~ 2.2 



90, so that Xf /g if (xf) 1 ^ 2 
(with variations of order 10 %). Then, from Eq. Q 

(o-ann v) > 7.3 x 10" 38 cm 2 /(^cDM^ 2 )max- Using the 
value (ficDM^jmax = 0.13, we obtain (cr ann v) > 5.6 x 
10~ 37 cm 2 . If a > l/(2x/)|6|, this implies 



B. 



< c ann u > > 3 x 10 37 cm 2 . 



Upper bound on £ 2 < CT an n« >o 



(8) 



Often it turns out that < a ann v >o may be orders 
of magnitude larger than the lower limit of Eq. JSJ. 
However, this fact does not automatically imply very 
large values for the antiproton flux, since the relevant 
quantity which enters in the antiproton flux is not sim- 
ply < cr ann w > , but instead £ 2 < cr ann « >o, through 
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the factor T. Indeed, the quantity £ 2 < a ann v >o co- 



incides with < cr a 



x v > when fl x h 2 > (ficDM^minj 



but it is proportional to < (7 ann i> >o / (cann v) when 
tt x h 2 < (Q C d 



mm 



then it has a maximum at Q. x h 2 



(^CDMfr 2 )min |13Ul4l ]. Let us call 77 the set of parameters 
of the particle-physics model which describes our generic 
cold relic. The property, that we have just discussed, 
states that the maximum of the quantity £ 2 < <7 ann v >o: 



(£ 2 < CannW >o)max =< C ann 7J > (,,=,,' (9) 

occurs when the model parameters rj have values 77', such 
that (il x h 2 ) v=n > = (S!cDMfr 2 )min, that is (using Eq. JjJ), 
when 
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Using the estimate already employed above, i.e. 
Xf/g i ,(xf) 1 ^ 2 ~ 2.2, from Eq. I|10[) one obtains 



(< cr ann v >)r,=n' = 



(Q x h 2 ) t 



7.7 x 10" 37 cm 2 , 



(11) 

where, in the last step, the value (fi x /i 2 ) m i n = 0.095 
is used. Thus, within a factor of 2, the maximum of 
£ 2 < c ann v > is stable for all sets of parameters 77 which 
satisfy Eq. i|llf> and is given by 



(£ 2 < Cannf >o) 1: 



x 10 



" 37 cm 2 . 



(12) 



If we take the mass m x to be one of the model parame- 
ters, this limit is independent of m x , provided the other 
parameters vary within ranges which allow solutions of 
Eq. (jnj. 

As a consequence of the previous properties, we find 
that the maximum of the factor T, as a function of m x , 
is expected to decrease simply as T max oc ni~ 2 . 

We note that the above properties apply also to any 
other primary flux of cosmic particles due to WIMP self- 
interactions in the halo, e. g. to the gamma-ray or 
positron flux. 



III. PRIMARY COSMIC ANTIPROTON FLUX 
FROM SELF-ANNIHILATION OF LIGHT 
NEUTRALINOS 

Now we finalize our previous considerations to the case 
of light relic ncutralinos. 



A. The Supersymmetric Model 

The supersymmetric scheme employed here is an ef- 
fective Minimal Supersymmetric extension of the Stan- 
dard Model (MSSM) at the electroweak scale, where 
gaugino-mass unification is not assumed. This is de- 
fined in terms of a minimal number of parameters, only 
those necessary to shape the essentials of the theoreti- 
cal structure of MSSM and of its particle content. The 
assumptions that we impose at the electroweak scale 
are: a) all squark soft-mass parameters are degenerate: 
?7iq ; = rriq, b) all slepton soft-mass parameters are degen- 
erate: mj. = mf c) all trilinear parameters vanish except 
those of the third family, which are defined in terms of 
a common dimensionless parameter A: A~ b = A^ = Aniq 
and A* = Amj. As a consequence, the supersymmetric 
parameter space consists of the following independent pa- 
rameters: M2, ^, tan /?, mAi Triq, mj, A and R = M1/M2. 
In the previous list of parameters we have denoted by jj, 
the Higgs mixing mass parameter, by tan/? the ratio of 
the two Higgs v.e.v.'s, by rriA the mass of the CP-odd 
neutral Higgs boson, and by Mi, Ma the U(l), SU(2) 
gaugino masses, respectively. 

In the numerical random scanning of the supersym- 
metric parameter space we use the following ranges: 
1 < tan/3 < 50, 100 GcV < M 2 < 1000 GeV, 
100 GeV < m^mj < 1000 GeV, sign(/j) = -1,1, 
90 GeV < m A < 1000 GeV, -3 < A < 3, in addi- 
tion to the above mentioned range 0.01 < R < 0.5. We 
impose the experimental constraints: accelerators data 
on supersymmetric and Higgs boson searches, measure- 
ments of the b — ■> s + 7 decay and of the branching ratio 
B s — > /i + /_t~, and measurements of the muon anomalous 
magnetic moment a M = (g^ — 2)/2. For the ranges used 
for these observable and other details of the model we 
refer to Ref. 0. 



B. The supersymmetric flux factor T 

Figs. 2] HI El gi ye the scatter plots of the quantities 
< o~axm v >o> £ 2 < CT ann ?j >o and T versus m x , limited to 
supersymmetric configurations which satisfy the approx- 
imation a > l/(2xf)\b\ employed above. (Red) crosses 
denote configurations with a relic abundance in the cos- 
mologically relevant range 0.095 < fl x h 2 < 0.13, (blue) 
dots denote susy configurations where rescaling is effec- 
tive (i.e. neutralinos form a subdominant species of relic 
particles) . In Fig. ^ we notice the effect of the lower 
bound on < o~ antl v >q, implied by the cosmological up- 
per bound on Ocdm^ 2 (see Eq. iJSJ). The rapidly ris- 
ing of the scatter plot as m x reaches the value of 45 
GeV is due to the self-annihilation process through the 
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FIG. 1: Scatter plot of < (7 ann u >o vs m x . Red crosses denote FIG. 2: Scatter plot of £ 2 < a ann v >o vs m x . Notations as in 
the supersymmetric configurations whose relic abundance is Fig. 
in the range 0.095 < fl x h 2 < 0.13, while blue dots denote 
configurations with Q, x h 2 < 0.095. 



Z-boson exchange, superimposed to a similar enhance- 
ment extending to higher values of m x , originating in the 
self-annihilation process with the exchange of the light- 
est CP-even neutral Higgs boson h. The upper frontier 
of the scatter plot at low values of m x is determined by 
the experimental lower bound on the mass of this Higgs 
boson. 

Fig. [21 displays the upper bound on £ 2 < a ann v >o, 
whose approximate estimate is given in Eq. H12|) . 

Finally, we note that the upper frontier of the scatter 
plot for T in Fig. [3] clearly displays the simple behavior 
TTmax oc rn~ 2 , as derived in Sect. Ill Bl The lower part 
of the plot is composed by configurations with a large 
rescaling in the local density. 

C. The differential antiproton spectrum dNp/dTp 

To evaluate dNp/dTp, we follow the treatment of Ref. 
||. In case of neutralino masses below the thresholds 
for gauge-bosons, Higgs-bosons and t quark production, 
antiprotons originate from the hadronization into p's of 
the quark and gluon pairs produced in the neutralino self- 
annihilation. For light ncutralinos (m x < 50 GeV), which 
are mainly binos with a slight mixing with a higgsino 
component the dominant final states are the ones 
into bb and into t~t + , the channel into bb being largely 
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FIG. 3: Scatter plot of the quantity T vs m x . Notations as 
in Fig. □ 

prominent for m x < 25 GeV. This property is displayed 
in Fig. Our calculation of the energy spectra has been 
performed by using a Monte Carlo simulation with the 
PYTHIA package fl5| . 

For neutralino masses which kinematically allow other 
final states (gauge-bosons, Higgs-bosons and ti pairs), 
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FIG. 4: Scatter plots of the branching ratios for the neutralino self-annihilation into bb (left panel) and r + r vs m x (right 
panel) . 



the full decay chain down to the production of quarks 
and gluons has been evaluated analytically. The final 
antiproton spectrum is then calculated from the previous 
results by boosting the differential energy distribution to 
the rest frame of the annihilating neutralinos. Details of 
our procedure are given in Refs. 0, 111] ]. A sample of 
results is shown in Fig. [5] where the spectra of antipro- 
tons per annihilation event are shown for a sample of the 
neutralino masses and self-annihilation into bb. 

D. Cosmic ray propagation in the Galaxy 

As they propagate throughout the Galaxy, charged 
cosmic rays mostly bounce on the spatial irregularities 
of its magnetic fields. That process is well described by 
mere diffusion with the energy dependent coefficient 

K = K PK S , (13) 

where 1Z stands for the particle rigidity. The acceleration 
of primary species by supernovae driven shock waves as 
well as their subsequent interactions with the interstel- 
lar gas take place in a thin galactic disk that is sand- 
wiched above and beneath by two large diffusion layers 
with thickness L. Because the magnetic field inhomo- 
geneities - actually the Alfven waves - move with speed 
Va, cosmic rays undergo inside the disk some diffusive 
reacceleration that come into play with the ionization, 
Coulomb and adiabatic energy losses. Particles are fi- 
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FIG. 5: The differential antiproton spectrum dNp/dTp for 
annihilation in the bb channel, as a function of the antiproton 
kinetic energy Tp. From left to right: ?n x = 10, 20, 30, 40 and 
100 GeV. 



nally wiped away by a vertical convective wind whose ve- 
locity is V c . The propagation of charged cosmic rays can 
be well accommodated in a cylindrical two-zone diffusion 
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model. The particle radial abundances may be expanded 
as series of Bessel functions Jo(aj r/R ga {) where a* is the 
i-th zero of the function Jo and where i? ga i = 20 kpc is 
the radius of the propagation region. For a complete de- 
scription of the semi-analytic code on which the present 
analysis is based, we refer the reader to Ref. |l6 U17U18l| . 
Suffice it to say that the five cosmic ray propagation pa- 
rameters mentioned above - namely the diffusion coef- 
ficient normalization Kq and index (5, the confinement 
layers thickness L and the velocities V c and Va ~ have 
been constrained by comparing the flux predictions for 
various cosmic ray species with observations. The most 
stringent observable is the boron to carbon ratio B/C - a 
typical secondary to primary relative abundance - whose 
analysis within our diffusion model has been detailed in 
Ref. [IRIl7j. The relevance of the latter has been further 
established by the compatibility of the B/C results with 
several other observed cosmic ray abundances [lTlHiLlioj . 

As a matter of fact, the conventional background for 
antiprotons is produced in the galactic disk by the spalla- 
tion of proton and helium cosmic rays on the interstellar 
medium. A calculation of the flux at the Earth of that 
population of secondary antiprotons has been performed 
in Ref. [lflj with the same propagation model as for the 
B/C analysis. The secondary flux has been derived con- 
sistently by employing the propagation parameters that 
have been obtained from the full and systematic anal- 
ysis of Ref. 0| on stable nuclei. The interstellar (IS) 
fluxes have been solar-modulated according to the force 
field approximation in order to obtain top-of-atmospherc 
(TOA) fluxes. Throughout this paper, we compare our 
results with data taken during minimal solar activity pe- 
riods and we fix the modulation potential <j> to 500 MV. 
It was found that the agreement between low-energy an- 
tiproton data and the estimate is excellent. Furthermore, 
the band within which the secondary flux lies is fairly 
restrained as shown in Fig. [5] In Ref. it was actu- 
ally inferred a ~ 20% astrophysical uncertainty - related 
to the propagation model - of the same order of magni- 
tude as the theoretical error arising from nuclear physics. 
Taking adiabatic energy losses into account results - be- 
low a kinetic energy of 1 GeV - into a small increase of 
our new secondary flux with respect to the previously 
published spectrum of Ref. 0] . 

Constraining cosmic ray propagation with B/C pro- 
vides a strong handle on any species that originates from 
the disk of the Milky Way and allows in particular the as- 
trophysical uncertainties on secondary antiprotons to be 
mild as discussed above. At variance with that secondary 
component, primary antiprotons are produced by the an- 
nihilations of neutralinos spread all over the galactic halo. 
As extensively discussed in Ref. 0], the corresponding 
uncertainties span now two orders of magnitude, when 
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FIG. 6: The secondary antiproton flux lies within the uncer- 
tainty band - delineated by the two solid black lines - which 
corresponds to all the possible propagation schemes that are 
compatible with the B/C data. Experimental observations at 
solar minimum are featured, in good agreement with the the- 
oretical estimate. Full circles: BESS 1995-97 |2lll : open squares 
BESS 1998 0; stars: AMS [3; empty circles: caprice E3| . 



the propagation parameters are varied between the min- 
imal and maximal configurations presented in Tab. 
Primary fluxes depend sensitively on the thickness of the 



case 


8 


Ko 


L 


Vc 


Va 






[kpc 2 /Myr] 


[kpc] 


[km s" 1 ] 


[km s _1 ] 


max 


0.46 


0.0765 


15 


5 


117.6 


med 


0.70 


0.0112 


4 


12 


52.9 


min 


0.85 


0.0016 


1 


13.5 


22.4 



TABLE I: Astrophysical parameters compatible with B/C 
analysis and yielding the maximal, median and minimal pri- 
mary antiproton flux. 



confinement layers and also on the convective wind that 
wipes cosmic ray species away from the galactic disk. It 
is not surprising therefore if the largest primary antipro- 
ton yield corresponds to the combination L = 15 kpc and 
V c = 5 km s _1 , whereas the smallest flux obtains when 
L = 1 kpc and V c = 13.5 km s _1 . The Alfvenic velocity 
Va is strongly correlated with the normalization constant 
K because the B/C ratio determines the diffusive reac- 
celeration parameter Kee oc V^/K - see Ref. - 
with an accuracy of ±15%. The parameter Va will not 
be further discussed. Finally, we remind that - as ob- 
tained in Refs. 0, [n| - the specific form assumed for the 
dark matter distribution in the Galaxy is fairly irrelevant 
in the calculation of the propagated primary antiproton 
flux. In particular, for not too thick confinement lay- 
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ers and strong convection winds, the cosmic ray diffusion 
range is small and solar circle abundances are blind to 
the center of the Milky Way and its putative neutralino 
cusp. 

IV. RESULTS AND CONCLUSIONS 

Because our analysis is focused on light neutralinos - 
m x < 40 GcV - our main interest is in low-energy an- 
tiprotons. That is why we have selected the observations 
performed in the GeV re gion by Bess 1995-97 [2l|, Bess 
1998 [13] and AMS 1998 l^. All these experiments op- 
erated at solar minimum for which <f> = 500 MV. The 
corresponding 32 data points are presented in Fig. Eland 
are already very well explained by a pure secondary com- 
ponent. Actually, we have derived a x 2 °f 33.6 in the 
case of the median configuration of Tab. [I] As discussed 
in Sec. IIII Dl should the galactic cosmic ray propagation 
parameters be varied so as to keep the B/C reduced x 2 
below a generous value of 1.8, the secondary antiproton 
flux would change by less than 20% - see Ref. [l9j. To 
put our discussion on a quantitative basis, we need to es- 
tablish a criterion for deciding whether or not a specific 
supcrsymmctric configuration is excluded by the above 
mentioned antiproton data. To this purpose, we compute 
the primary flux which such a configuration yields, add it 
to the secondary component and derive the correspond- 
ing x 2 - Notice that Caprice has been disregarded 
in our calculation as it operated at a slightly higher en- 
ergy than the range in which we are interested. Then 
we decide that a specific supersymmetric configuration is 
excluded if the corresponding x 2 exceeds a critical value 
of 60. For 32 degrees of freedom, this translates into a 
disagreement at the 99.5 % C.L. It is worth stressing that 
our exclusion criterion is purposely taken on the conser- 
vative side. 

A final word of caution is in order before we proceed. 
The problem is complicated by the fact that the galactic 
cosmic ray propagation model is not unique. Actually, in 
the present analysis, we have varied the five propagation 
parameters presented in Sec. IIII Dl and considered all the 
possible combinations that are in good agreement with 
stable nuclei 0|. We keep the conservative attitude of 
selecting propagation models as long as they lead to a 
reduced x 2 e dB/c — ^ on the B/C data. 

To commence, we focus our investigation on the su- 
persymmetric configurations which the present data on 
low-energy cosmic ray antiprotons constrain most. Those 
configurations are associated to large antiproton yields 
and correspond to high values of the effective annihila- 
tion cross-section £ 2 < a ann v > - We concentrate here on 
cases A) and B) for which the neutralino relic density is 



relevant to cosmology. Case A) corresponds to the largest 
antiproton signal and consequently to the strongest con- 
straint as is clear in the panel a of Fig. [7| where the x 2 
is plotted as a function of the diffusive halo thickness L. 
(Black) squares, (green) big circles, (red) smaller circles 
and (blue) dots correspond to a neutralino mass m x of 
10, 20, 30 and 40 GeV, respectively. The constellation of 
points of the same shape (and color) - which corresponds 
to a specific neutralino mass - is obtained by scanning 
the entire set of cosmic ray propagation models that are 
compatible with B/C - as discussed above. The effective 
annihilation cross-section £ 2 < tr ann w >o reaches here a 
maximal value of 10 -36 cm 2 as may be readily inferred 
from the scatter plot of Fig. |21 and its upper boundary. 
Case B) corresponds to the minimal value of the anni- 
hilation cross-section < a ann v >o that is yet compatible 
with WMAP. The neutralino relic abundance is set to 
the upper bound (ricDM)max^ 2 = 0.13. The correspond- 
ing scatter plot is presented in panel b of Fig. [7| The 
lower horizontal line at x 2 = 33.6 corresponds to the fit 
with the secondary antiproton component alone. The up- 
per horizontal line indicates the discriminating x 2 value 
above which the fit to the low-energy antiproton data is 
no longer acceptable. The large spread exhibited by the 
scatter plots of Fig. [3 illustrates the strong sensitivity of 
the primary antiproton flux to the astrophysical parame- 
ters that account for cosmic ray propagation throughout 
the Milky Way and its halo. This is at variance with 
the stability of the secondary flux as already remarked 
above. As is clear from panels a and b, light neutrali- 
nos of cosmological interest - i. e. with a relic abundance 
in the WMAP range - are compatible with the present 
data on cosmic ray antiprotons only if the parameter L is 
on the very low side of its physical interval. Should the 
diffusive halo thickness L turn out to be larger than 2 
kpc - a quite realistic assumption given the existence of 
electron synchrotron radiation far above the galactic disk 

- neutralinos lighter than 20 GeV would be excluded. If 
we now assume that L > 5 kpc, we readily conclude that 
neutralinos with sizeable relic abundance must be heavier 
than 40 GeV. In order to illustrate how our x 2 exclusion 
criterion translates on the observations, we have plotted 
in Fig.|Hlthe total antiproton spectrum - (red) solid curve 

- for a 20 GeV neutralino with maximal relic abundance. 
The median configuration of Tab. |U has been selected for 
the cosmic ray propagation parameters. It corresponds 
to the best fit to the B/C data. The primary (blue long- 
dashed curve) and secondary (blue dotted line) compo- 
nents add up to yield an antiproton signal well in excess 
of the observations. The corresponding x 2 reaches actu- 
ally the unacceptable value of 200. 

The last panel of Fig. features supersymmetric con- 
figurations for which the effective annihilation cross- 
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FIG. 7: The \ 2 is featured as a function of the diffusive halo thickness L. Black squares, green big circles, red smaller circles 
and blue dots correspond to a neutralino mass m x of 10, 20, 30 and 40 GeV, respectively. At fixed neutralino mass m x , 
each point is associated to a specific combination of the galactic propagation parameters yielding a good fit for the B/C ratio. 
Panels a, b and c respectively correspond to cases A), B) and C) defined in the text. The upper horizontal line indicates the 
discriminating \ 2 value above which the fit to the low-energy antiproton data is no longer acceptable. The lower horizontal 
line corresponds to the \ 2 calculated with the secondary component alone for the median configuration of Tab. [I] 



section £ 2 < cr ann v >o is minimal and rescaling is max- 
imal. The antiproton signal is at its weakest level and 
we do not expect low-energy antiproton data to be very 
constraining. Panel c of Fig.0indicates nevertheless that 
a 10 GeV neutralino is excluded provided that the diffu- 
sive halo thickness L exceeds 2.5 kpc. Should L be larger 
than 10 kpc, the limit rises to 20 GeV. That panel shows 
also to what extent cosmologically subdominant neutrali- 
nos with a mass in the range 20 GeV < m x < 40 GeV 
escape conflict with present data. A potential distor- 
tion is difficult to unravel from the antiproton spectrum 
given the available observations because error bars are 
too large. Actually, Fig. [5] and ^| feature respectively 
the primary antiproton flux of a 10 and 20 GeV neu- 
tralino (blue long-dashed curves) in the case C) of max- 
imal rescaling together with the conventional secondary 
component (blue dotted lines). The cosmic ray propa- 



gation configurations that have been extracted from the 
panel c of Fig.0- a (black) square if the neutralino mass 
is 10 GeV or a (green) big circle when m x is twice as 
large - lie on the upper horizontal line for which the 
fit to the low-energy antiproton data yields a \ 2 value 
of 60. Below an antiproton kinetic energy of ~ 1 GeV, 
the global signal represented by the (red) solid curves is 
slightly shifted upwards with respect to the secondary 
spectrum whereas above 1 GeV, both spectra are identi- 
cal. Richer samples of experimental data are a necessary 
condition for disentangling a supersymmetric signal from 
background secondary antiprotons. The spectrum of the 
latter suffers nevertheless from theoretical uncertainties - 
cosmic ray propagation and nuclear cross-sections - that 
need to be further reduced in order to make antiproton 
measurements a useful probe for cosmologically subdom- 
inant ncutralinos in the 10 to 30 GeV range. 
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FIG. 8: The primary antiproton spectrum - blue long-dashed 
curve - has been computed for a 20 GeV neutralino and corre- 
sponds to the case B) of a maximal relic abundance. The me- 
dian configuration of Tab.UJhas been assumed. When added 
to the secondary component - blue dotted line - the primary 
signal results into a total antiproton flux - red solid spectrum 
- well above the data points. We actually infer a x 2 value of 
200 to be compared to the exclusion limit of 60. Symbols for 
experimental data as in Fig. 
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FIG. 9: The antiproton spectrum is featured together with 
its primary - blue long-dashed curve - and secondary - blue 
dotted line - components for a 10 GeV neutralino. The super- 
symmetric configuration corresponds to maximal rescaling - 
case C) - and the specific set of astrophysical parameters that 
has been extracted from panel c of Fig. Q This case corre- 
sponds to x 2 =60, and then is marginally acceptable. Symbols 
for experimental data as in Fig. |S| 

It is very instructive to pursue the analysis of the max- 
imal rescaling case C) in terms of other astrophysical pa- 
rameters. The three scatter plots of Fig. ^2 feature the 
same selection of astrophysical configurations drawn from 
Fig.0- panel c - for which the fit to the antiproton data 
is acceptable with \ 2 < 60. Big (black) circles, smaller 
(green) circles and (red) dots correspond to a neutralino 
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FIG. 10: Same as in Fig^with a 20 GeV neutralino. 

mass m x of 10, 20 and 30 GeV respectively. The smaller 
the latter, the more abundant the neutralinos at fixed 
mass density and the stronger the antiproton annihila- 
tion signal. This trend is clear in each of the panels of 
Fig. It is even clearer in panel c. The relevant annihi- 
lation cross-section is in that case the minimal value of 
£ 2 < c ann u >o which increases rapidly with decreasing 
neutralino mass - see the lower boundary of the scat- 
ter plot in Fig. |3 Low m x configurations are allowed if 
the primary antiproton excess which they yield does not 
propagate efficiently from the halo of the Milky Way to 
the solar system. Thin diffusive halos and strong galac- 
tic winds are preferred. This is particularly obvious for 
a 10 GeV neutralino. Notice how the big (black) circles 
concentrate in the upper-left corner of the V c versus L 
diagram. Remark also how the normalization Kq of the 
diffusion coefficient scales with the diffusive halo thick- 
ness L. In the second panel of Fig. El the ratio Kq/L lies 
actually in the range between 10~ 3 and 10 -2 kpc Myr -1 
irrespective of L. That ratio turns out to be crucial in 
the diffusion equation. Finally, we do not find any par- 
ticular correlation between 5 and L even if small values 
for the latter are still pointed towards in the last panel 
of Fig. EU 

To conclude our investigation, we would like to assess 
the potential of low-energy antiproton measurements for 
discovering heavier neutralinos. For that purpose, we 
present in the \ 2 versus L plane of Fig. ^| the same 
kind of scatter plot as in panel a of Fig. [7] Neutralino 
masses of 100 GeV (red stars) and 200 GeV (blue crosses) 
have been displayed. In the supersymmetric configura- 
tion which we have selected for each neutralino mass, the 
effective annihilation cross-section £ 2 < <J ann v >o reaches 
its maximal value - case A) - so that the antiproton sig- 
nal is the strongest. At fixed m x , each star or cross is 
associated to a combination of galactic propagation pa- 
rameters yielding a good fit for the B/C ratio. Notice how 
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FIG. 11: The three panels feature the same selection of astrophysical configurations from panel c of Fig.|7|for which the fit to 
the antiproton data is acceptable with x 2 < 60. Big black circles, smaller green circles and red dots respectively correspond to 
a neutralino mass m x of 10, 20 and 30 GeV. The panels successively display the cosmic ray propagation parameters V c , Ko/L 
and 5 as a function of the diffusive halo thickness L. 
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FIG. 12: Same as in Fig.Q panel a) but for m x =100 GeV - 
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FIG. 13: The % 2 is plotted as a function of the neutralino 
mass. Case A) has been selected with astrophysical param- 
eters leading to the maximal signal. Horizontal lines are the 
same as in Fig. [7| 



the constellation of (blue) crosses lies between the two 
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horizontal lines. Wc readily conclude that whatever the 
cosmic ray propagation model is, low-energy antiproton 
data are useless as long as m x > 200 GeV. Our statement 
is supported by Fig. 1131 where the maximal possible value 
of x 2 nas been plotted as a function of neutralino mass. 
The effective annihilation cross-section £ 2 < u ann v > 
has been once again tuned to its maximal value - case 
A) - and the astrophysical configuration max of Tab. [I] 
has been chosen. The x curve drops below the exclusion 
value of 60 when the neutralino mass exceeds 200 GeV. 
Observations at significantly higher energies will proba- 
bly be necessary to explore the regime where neutrali- 
nos are heavy. We leave such an investigation for a fu- 
ture publication but we cannot resist noticing that the 
Caprice data above 10 GeV already exhibit an excess 
with respect to the secondary antiproton background. 
If such a distortion is confirmed - by the forthcoming 
Pamela satellite mission j2j| for instance - it should have 
to be explained ^lj ■ Pamela will actually probe the an- 
tiproton spectrum between 80 MeV and 190 GeV and 
collect data during three years starting at the end of 
2005. As well as the Pamela experiment, other space 
missions are planned - such as AMS-02 on board the In- 
ternational Spatial Station 1261 o r the balloon-borne Bess 
Polar mission in Antarctica |2jj - or in project - such as 
the GAPS satellite experiment |2^|. They will measure 
the low-energy antiproton spectrum with very high ac- 
curacy, and their results will be of great relevance for 
improving the study developed in the present paper. Im- 
provements in the understanding of the propagation of 



galactic cosmic rays are foreseeable in the next future, 
thanks to long and ultra-long duration balloon missions 
[29l | and space-based expriments [2q|. Even if the best 
expectations are deserved to the determination of the 
diffusion coefficient power spectrum 5, data on the B/C 
quantity in the GeV/n region will also constrain the diffu- 
sive halo thickness L which, among the astrophysical 
parameters, plays the most important role in our previ- 
ous analysis. 

Finally, we recall that the conclusions of the present 
paper are drawn for the case of a smooth distribution of 
dark matter in the galactic halo. Should the dark matter 
distribution have some dumpiness, the antiproton signal 
would be enhanced j^H, and consequently some of the 
constraints on the previously discussed supersymmetric 
configurations would become more stringent. 
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